na sa cf\- | (j,-) 95^ 


NASA Contractor Report 167959 NASA-CR-1 67959 

19830010400 


USER'S MANUAL FOR THREE-DIMENSIONAL ANALYSIS 
OF PROPELLER FLOW FIELDS 


Denny S. Chaussee and Paul Kutler 

Flow Simulations, Inc. 

Sunnyvale, California 


January 1983 


i-i 


Prepared for 



lancu:, r,; : 
L' t 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Lewis Research Center 
Under Contract NAS 3-22375 



i 


TABLE OF CONTENTS 

Section Page 


1. SUMMARY 1 

2. INTRODUCTION 2 

3. LIST OF SYMBOLS 6 

4. ANALYSIS 10 

4.1. Derivation of Governing Equations 10 

4.2. Numerical Method 17 

4.3. Geometry and Mesh Generation Descriptions. . 23 

4.4. Boundary and Initial Conditions 36 

4.4.1. Tangency Condition 37 

4.4.2. Kutta Condition 41 

4.4.3. Subsonic Outflow 45 

4.4.4. Initial Conditions 54 

4.5. Output Data Manipulation 56 

4.5.1. Blade Properties 62 

4.5.2. Nacelle Properties 72 

5. MESH GENERATION PROGRAM 73 

5.1. Program Organization and Interaction .... 73 

5.2. Flow Chart 75 

5.3. Program Input 76 

5.3.1. Dictionary of Input Variables. ... 76 

5.3.2. Input Data Format 79 

t 

-12 

n- l 5 


T 


11 

Section Page 

6. THREE-DIMENSIONAL PROP-FAN FLOW FIELD PROGRAM. . . 83 

6.1. Program Organization and Interaction .... 83 

6.2. Flow Chart 84 

6.3. Program Input 85 

6.3.1. Dictionary of Input Variables. ... 85 

6.3.2. Input Data Format 87 

6.4. Program Output 88 

6.5. Operating Procedures and Trouble Shooting 

Hints 90 

7. DATA REDUCTION PROGRAM 91 

7.1. Program Organization and Interaction 91 

7.2. Flow Chart 92 

7.3. Program Input 93 

7.3.1. Dictionary of Input Variables .... 93 

7.3.2. Input Data Format 94 

7.4. Program Output 95 

7.4.1. Output Format 96 

7.5. Operating Procedures and Trouble Shooting 

Hints 97 

8. PROGRAM SUMMARIES 98 

8.1. Machine Compatibility 98 

8.2. Run Time Estimates 99 

REFERENCES 101 



Ill 


Section Page 

APPENDICES 

A. Sample Input 103 

A.l. Mesh Generation Program 103 

A. 2. Three-Dimensional Flow Field Program .... 103 

A. 3. Data Reduction Program 103 

B. Sample Output 115 

B. l. Mesh Generation Program 115 

B.2. Three-Dimensional Flow Field Program .... 115 

B.3. Data Reduction Program 115 





USER'S MANUAL FOR THREE-DIMENSIONAL 


ANALYSIS OF PROPELLER FLOW FIELDS 

by 

Denny S. Chaus see 
Paul Kutler 

Flow Simulations, Inc. 

1. SUMMARY 

A detailed operating manual is presented for the prop-fan 
computer code (in addition to supporting programs) recently 
developed by Kutler, Chaussee, Sorenson, and Pulliam while at 
the NASA's Ames Research Center. This code solves the 
inviscid Euler equations using an implicit numerical procedure 
developed by Beam and Warming of Ames. A description of the 
underlying theory, numerical techniques, and boundary conditions 
with equations, formulas, and methods for the mesh generation 
program (MGP) , three-dimensional prop-fan flow field program 
(3DPFP) , and data reduction program (DRP) is provided, 
together with complete operating instructions. In addition, 
a programmer's manual is also provided to assist the user 
interested in modifying the codes. Included in the programmer's 
manual for each program is a description of the input and 
output variables, flow charts, program listings, sample input 
and output data, and operating hints . 
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2 . INTRODUCTION 

A variety of advanced propellers is being designed for 
high-speed aircraft of the future. It has been shown from 
experimental tests that the new propeller termed the prop-fan 
is less noisy and more efficient than other models tested, 
and could result in saving a large amount of airplane fuel. 

The prop-fan is a small-diameter, highly loaded, multibladed, 
variable pitch propulsor (see Fig. 1) . The blades, which are 
twisted and tapered, utilize thin airfoil sections with tip 
sweep. They are integrated with a spinner and nacelle shaped 
to reduce the axial Mach number through the blading. Interest 
in the prop-fan propulsion system concept arose as a result of 
its potential for significant fuel savings over conventional 
high bypass turbofans . Interest from a computational point 
of view was stimulated by the need for understanding the flow 
phenomena around propellers and the potential for a more 
efficient computer generated design. 

Over the years, computational procedures for solving fluid 
flow problems provided an inexpensive but accurate means of 
determining the aerodynamic characteristics of complex 
configurations. In addition, they have provided the designer 
with an effective tool for maximizing aerodynamic efficiency 
without the expense of actually building and testing numerous 
designs. Finally, computational methods have often been 
capable of providing information not readily obtainable from 
experiments . 


Figure 1. - Prop-fan mounted in wind tunnel. 
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The prop-fan is a good example of the type of configuration 
for which it is difficult to experimentally obtain aerodynamic 
information needed for performance analyses and design. The 
blades are virtually impossible to adequately instrument 
because of the high rotational velocities, and their relatively 
small thickness. Also, details of the surrounding flow field 
can only be obtained by performing costly and time consuming 
flow field surveys. Computationally, however, the entire flow 
field can be determined from a single solution of the governing 
equations including near- and far-field effects which can be 
used in acoustic analysis programs and blade surface pressure 
distributions which can be used in structural and aerodynamic 
design analysis programs. 

The theoretical development and subsequent numerical 
solutions were concerned with simulating the inviscid flow 
about the prop-fan model used in the early experiments, 
i.e., the standard eight-bladed prop and spinner mated with 
an axisymmetric nacelle or afterbody (instead of the conventional 
three-dimensional' nacelle with inlet). Such a configuration 
thus requires that only the flow about a single blade be 
computed because of periodicity. The new and unique feature 
of the present work over conventional efforts is that the 
blades are unshrouded, i.e., the flow is permitted to spill 
over the blade tip. In effect, the entire flow field about 
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a lifting wing-body combination with the added complexity of 
rotation was computed. 

This report is broken down into two basic categories. 

The first is the Analysis section in which the theory used 
to develop the prop-fan flow field program is explained. This 
includes the derivation of the governing equations, a description 
of the numerical algorithm, body definition and mesh generation, 
boundary and initial conditions, and data reduction theory. 

The second part of this report describes the three computer 
programs that were developed to analyse the prop -fan flow field. 
These include the mesh generation program (MGP) , three- 
dimensional prop-fan flow field program (3DPFP) , and data 
reduction program (DRP) . The description of each program 
includes its organization and interaction with the other 
programs, a flow chart, sample input and output data, operating 
procedures and trouble shooting hints, and a basic description 
of each main program and subroutine. These computer programs 
were written in FORTRAN IV source language and developed for 
use on a CDC 7600 computer using the OPT = 2 compiler. 
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3. LIST OF SYMBOLS 
defined in Eq. (25) 

axial force or Jacobian matrix, 3E/3Q 

reference area, nr^ 

b 

blade width or chord 
defined in Eq. (25) 

Jacobian matrix 3F/9Q 
speed of sound, = (yp/p)^ 
defined in Eq. (25) 

Jacobian matrix 3G/9Q 
axial force coefficient 
force coefficient 
moment coefficient 
normal force coefficient 
pitching moment coefficient 

pressure coefficient also (Section 7.1) 

torque coefficient 

rolling moment coefficient 

side force coefficient 

thrust coefficient 

power coefficient 

yawing moment coefficient 

Courant number (see Eq. (7)) 

diameter of propeller 

drag, Eq. (95) 
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LIST OF SYMBOLS (Continued) 


dF 

dM 

dN 

dS 

e 

E 

F 

g 

G 

H 


H 


H. 


(iz.i-r,^) 


(i,j ,k) 

J 

L 

LEA 

S 

M 




M 


R 


Hy 


incremental force vector, Eq. (72) 
incremental pitching moment 
incremental normal force, Eq. (76) 
incremental side force, Eq. (76) 
total energy per unit volume 
vector of flux quantities in ^-direction 
vector of flux quantities in q -direction 
integrand of Eq. (57) 

vector of flux quantities in ^-direction 

vector of source terms created by cylindrical 
generalized coordinate transformation 

rothalpy 

total enthalpy 

unit normal vector of cylindrical coordinate 
system 

unit vectors of cartesian coordinate system 

Jacobian or Advance ratio (u ro /nD) 

lift 

leading edge alignment 
reference length 
Mach number 
pitching moment 
rolling moment 
yawing moment 
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n 

A 

n 

A 

n c 

A 

n f 

n l/2 

N 

P 

k ± , i=0 , 3 

^00 

q i ,i=l,5 
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-> 
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*2 

r 3 

R 4 
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s 

s^. 

(t,z,r,<j>) 
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LIST OF SYMBOLS (Continued) 

rotational speed 

unit normal, Eqs.(32) and (54) 

Eq. (33) 

Eq. (34) 

Eq. (38) 
normal force 
static pressure 
generalized metrics 
dynamic pressure = O.Sp^v 2 
components of Q 

torque or vector dependent variable of integration 

distance from coordinate system origin to 
incremental area d£, Eq. (80) 

right-hand side of Eq. (54) 

right-hand side of Eq. (55) 

right-hand side of Eq. (56) 

right-hand side of Eq. (57) 

constant surface. Figure 9 

side force 

entropy Eq . (22) 

cubic function used in initial startup Eq. (71a) 
independent variables, cylindrical coordinates 
thrust 

used in initial startup Eq. (71a) 
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LIST OF SYMBOLS (Concluded) 


u 

U, u 

y 

V 

V, v 


w 

W, w 
(x,y,z) 

Y 

Y i ,i=l,5 

x cp 

AB 

Ap 

AC, An , Ac 
P 

(?,£,n ,0 

e _ 


w 

Subscripts 

b 

i 

j 


physical velocity in z-direction 
contravariant velocity defined in Eq. (2) 
see Eq. (5) 

physical velocity in y- or r-direction 
contravariant velocity defined in Eq. (2) 

AAA 

velocity vector ui + vj + wk 
See Eq. (5) 

physical velocity in x- or ^-direction 
contravariant velocity defined in Eq. (2) 
cartesian coordinate system 
ratio of specific heats 

eigenvalues of gas-dynamic equations (see Eq. (6)) 

center of pressure 

blade twist 

incremental pressure 

computational mesh spacing 

density 

transformed independent variables (Eq. (1)) 
explicit smoothing coefficient in implicit algorithm 
implicit smoothing coefficient in implicit algorithm 
angular velocity 

body 

integer mesh point location in ^-direction 

integer mesh point location in n-direction 
integer mesh point location in ^-direction 

free stream conditions 
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4. ANALYSIS 

4.1. Derivation of Governing Equations 

To enhance numerical accuracy and efficiency, non- 
orthogonal coordinate transformations of the governing 
equations in a particular base coordinate system are employed. 
This maps the surface of the nacelle and both sides of the 
blade onto constant coordinate surfaces which facilitates the 
application of boundary conditions and permits grid point 
clustering at the body when the dependent variables are 
expected to undergo rapid changes. Use of such transformations 
in addition permits utilization of uniform discretization 
formulas and well-ordered interior grid point solution 
algorithms. Under this transformation the equations can 
still be written in conservation law form to take advantage 
of the shock capturing properties. 

The basic orthogonal coordinate system utilized is 
cylindrical with z oriented along the rotational axis, 
r extending radially outward from the z-axis, and <f> the 
meridional angle measured from a vertical plane (see Fig. 2) . 

It should be reiterated that for this study, only the flow 
between two of the blades is computed, i.e. between the face 
side of one blade and the camber side of the other blade 
because of periodicity. The cylindrical coordinates are 



Figure 2. Coordinate System 
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transformed to align the blade and nacelle surfaces with 
various computational planes according to the following: 

T — t 

5 » £(t, z, r, cf>) 

n = n(t,z,r,cf>) (1) 

5 “ C(t,z,r,(J)) 

This generalized nonorthogonal coordinate transformation maps 
the spinner and nacelle onto a constant n -plane and each side 
of the blade, i.e., the camber and face sides, onto parts of 
a constant £ -plane. The remaining parts of the constant 
£ -planes are periodic surfaces. The radial far-stream and 
outflow boundaries are situated far enough from the prop-fan 
to minimize the reflection of waves. 

The governing partial differential equations in weak 
conservation law form for cylindrical coordinates under the 
assumptions of inviscid compressible flow and a perfect, 
non-heat conducting gas for the transformation given in 
Eq. (1) are: 


Q T + E e + F n + C ? + H = ° (2) 

where 


p 


PU 

pu 


pulH- P 5 z 

pv 

• e = j 

pvW-p£ r 

pw 


pwWp^/r 

e 


(e+p)U-p£ t 
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pV 


pW 


pv 

P uV+pn z 


puW+pc z 


puv 

puV+pn r 

• G = i 

pvW+pc^ 

■ H = k 

p(v 2 -w 2 ) 

pwV+pn^/r 


pwW+pt /r 
<P 


2pvw 

(e+p)V-pn t 


(e+p)W-p? fc 


(e+p) v 


and 

U = 5 t + u£ z + v r + wi; /r 

V = Hf + un + v + wn./r 
t z r cp 

W = c + uc + v + w£ . /r 

t z r <j) 

U, V, and W are the contravariant velocities written without 

metric normalization. J is the transformation Jacobian and 

is defined below. Use of the Euler equations in conservation 

law form guarantees the accurate calculation of the shock 

waves occurring near the blade tips. 

In the conservative variables of Eq. (2), the pressure 

p is nondimens ionalized by p^,, the density p by p^, and the 

cylindrical velocity components u, v, and w by a ra //y where 

2 

a ro is the free-stream speed of sound (a ra = Ypoo/Poo) and y is 
the ratio of specific heats. Other quantities made dimension- 
less are the time t by ta^/CD/y) and the angular velocity 
w by wD/y/a^. The pressure, density, and velocity components 


are related to the total energy per unit volume e by the 
following equation for an ideal gas : 

e = p/(y-l) + p(u 2 + v 2 + w 2 )/2 (3) 

The metrics of Eq. (2) are obtained by the chain rule 
expansion of z^ 

yield the following expressions: 


= 

Z t^z 

- Vr - 


^z - 

(r <j> 
v re 

- 4>_ r ) 

Y n e 

/X 

= 

" Z T n Z 

- r T n r - 


^r 

(4 Z 

e 

- z d> ) 

hV 

/I 

«t = 

' Vz 

‘ Vr * 


= 

<V 5 ■ 

- z r ) 

e n 

/I 

n z = 

<Ve - 

r^)/I 


? z = 

( Vn ' 

■ ♦sV 

/I 

n r = 

( Vs ' 



c r " 

( V? ' 

• V C > 

/I 

-e- 

II 

(r 5 Z C ' 

z r ) /I 

€ ? 


C <j> = 

(z 5 r n ' 

' r 5 z n> 

/I 


where 

I = 1/J = z r r d> + r d> z + * z r - z_d> r„ - r z d> - d>_r z 
i rn c ^ n c 5 n c £ ire n e 

The quantity 4> is the rotational or angular velocity of the 
propeller. 

In general, the metrics of Eq. (4) are not known 
analytically and must be determined numerically at each step 
of the intergration procedure. To accomplish this, second- 
order central-difference formulas are used at interior points 
and three-point one-sided formulas are used at the boundaries. 


'S’ 


etc. and solved for 


^ * Sc : 


etc. to 
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The implicit algorithm to be discussed below requires 
the analytic determination of the Jacobians A = 3E/3Q, 

B = 3F/3Q, and C = 3G/3Q. These Jacobian matrices resulted 
from the time linearization of E, F, and G in deriving the 
numerical algorithm. The Jacobian matrices can be written 
in general as follows : 


*0 


i 


"2 


r 


0 


I 


A, B. 
or C 


V-r q 7 

r 

j Ko + V + r^Z-T) ^ 
1 

1 

I - K (r -1) + ll !i | 

1 1 q l Y q x | 

*1 (Y-l) 



I 

1 

1 

1 

j K 0 + V + VZ-Y) ^ 
1 

1 1 

1 1 

I ^(Y-Z) + i 

1 2 q l r q l 1 

1 

*2 (Y-l) 

(5) 

r i 

1 

! i 

1 

1 

1 q 4 K 1 q Y 

K, — - =5(Y-1) -i 

1 2 q l r q l 

1 
1 

J 1 

| *0 + » + ^-Y> ^ I 

I 1 

K 3 

^ (Y-l) 


-T ^ V +<Y-1)V 2 

1 

| -<Y-l)jp 9 + KjY h 
j q l 1 q l 

1 

| -CY-i) 2a * + v 3i 

| !4 v+ _3 T i | 

k 0 + y V 



! - JjS v 2 

i -h?* 

! ' r <I r ) 72 ! 
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where V = (k^ + k 2 q 3 + q 4^ q l and ^ + + q 3 + q 4^/ q r 

The terms q^, i = 1, 5 are the components of the dependent 
variable of integration Q in Eq. (2). To obtain, for example, 

A, let k Q = £ t , k x = C 2 , k £ = C r and k 3 = ^ in Eq. (5) . 

As a measure of the productivity of the implicit algorithm 
to be described below compared to that of conventional explicit 
procedures, the value of the Courant number can be computed. 

It requires the eigenvalues of the matrices A, B, and C. 

They are as follows : 


a 


1,2,3 



+ uk^ + vk^ + f k 3 


and 

k 2 1/2 

a 4,5 = k 0 + ^1 + ^2 + r k 3 - C l k l + k 2 + ^ I 

where for A, k Q = C t » etc.; for B, k Q = n t » etc.; 

for C, k^ = £ t , etc.; and c equals the local speed of sound 

O 

(c^ = yp/p) • The Courant number is defined as follows: 


™ - 'WI'Wl < 7 > 

where it is assumed that AC = An = Ac = 1, and a max is the 
maximum value of the eigenvalues of all the nodal points . 

For the prop-fan problem described here, Eq. (2) is 
solved in a time-asymptotic fashion with interest only in 
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the steady-state solution. This is a result of the mesh 
rotating with the blade and the fact that a cylindrical 
coordinate system is used. The Q t term of Eq. (2) approaches 
zero as x becomes large, thus establishing a convergence 
criterion. 

4.2. Numerical Method 

The numerical algorithm used to solve the conservation- 
law form of the Euler equations is based on a class of 
completely implicit, noniterative, ADI (alternating direction 

implicit) schemes developed by Lindemuth and Killeen^, 

f\ o 4 5 

Briley and McDonald'* 1 * , and Beam and Warming ’ . The 

particular method is a generalization of a conservative, 

approximate factorization scheme in the "delta" form^ . 

The procedure has been successfully applied by Steger and 

6 7 8 

Kutler and Kutler' for inviscid flows and by Steger , 

a 

Pulliam and Steger , and numerous others, Refs. 10 and 13, 
for viscous flows. Use of the implicit procedure helps 
remove the stiffness of the problem introduced by a fine 
mesh. Thus, for this problem the implicit procedure permits 
an integration stepsize large enough to obtain steady-state 
solutions in considerably fewer iterations than conventional 
explicit procedures. 


T 
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As applied to Eq. (2) the implicit, spatially factored 
algorithm using Euler implicit time differencing takes the 
form 

1 1 + At<S^A n - e i (J" 1 V^J) n | 1 1 + A t<5^B n - J) n | 

1 1 + At6 C n - e (j'Vi J)“| (Q n+1 - Q n ) = 

^ 1 

-At 1 6 E n + 6 F 11 + H n + S G n | - ^ (J n ) I (V-A ) 2 

+ <V A ) 2 + (V A ) 2 | (JQ) n (8) 

? K n n 

where A, B, and C are the Jacobian matrices 9E/9U, 9F/9U, 

and 9G/9U, respectively, and I is the identity matrix. 

and 6 are second-order central-difference operators, and 

K n 

V and A represent the conventional forward and backward 

difference operators. The quantities multiplied by 

on the left-hand side of the equation represent implicit 

second-order smoothing terms (which were added later to the 
8 9 

basic algorithm ’ ) while the quantity multiplied by 
on the right-hand side represents an explicit fourth-order 
smoothing term. 

The solution of Eq . (8) consists of first forming the 
right-hand term (also called the steady-state or explicit 
part) at each grid point. Each of the implicit operators, 

(for instance {X + At 6^ A n - e^(J ^V^A^J) 11 } ) represent 
block- tridiagonal matrices which must be inverted sequentially 
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to obtain AQ n = Q n - Q n . A block lower-upper-decomposition 

12 

algorithm, see Isaacson and Keller , is used for the inversion 
process . 

For steady-state calculations AQ n = Q n "*"l - Q n 0 

and the solution then satisfies the steady-state finite- 

difference (right-hand side of Eq. (8)) equations. Linear 

stability analysis shows unconditional stability for the 

"delta" form of the implicit approximate factorization 

algorithm. In actual practice, though, time step limitations 

are encountered, although they are usually much less stringent 

than explicit stability bounds . The smoothing terms have been 

added to control nonlinear instabilities. Choices of the time 

step At and the smoothing coefficients (e ,e.) are usually 

dictated by experience. Linear analysis for the smoothing 

terms does show that for e. = 0, e is bounded (e < 1/4) 

i e e 

for stability, but for >> e , & is not constrained to 

any limit for stability. 

The metric terms are obtained using second-order finite- 
differences for terms such as, z^, in Eqs . (4). Fourth-order 
accuracy in the steady-state can be obtained in an efficient 
manner by introducing fourth-order finite differences for the 
convective terms on the right-hand side of Eq. (8), while 
retaining second-order differences on the left-hand side. For 
Euler implicit time differencing this is a stable and accurate 
process. For more details on the above, see Pulliam and 
Steger^ . 
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The major portion of the computational work in an implicit 

finite -difference algorithm is contained in the solution of 

the set of simultaneous equations, (i.e., the block tridiagonal 

equations) . When even an implicit algorithm is applied to a 

system of partial differential equations, one obtains block 

matrix -vector equations that are complicated and time-consuming 

to solve. A method for uncoupling the solution process through 

a diagonalization of the block -matrix structure has been 

13 

presented by Pulliam and Chaussee . The method was applied 
to an implicit approximate -factorization algorithm for the 
two- and three-dimensional inviscid Euler equations in general 
curvilinear coordinates . 

The Jacobian matrices. A, B, and C, have a set of 
eigenvalues and a complete distinct set of eigenvectors. 
Similarity transformations (see Warming, Beam, and Hyett 1 ^) 
can be used to diagonalize A, B, and C where 


A 



A = Ty.Aj.T~ 1 , B = T A T' 1 , 
IKK n n n * 

r % 

D [U.U.U.U + c(K 2 z + K 2 x + ^ /r z) , 

D [v , V , V , V + c(n‘ + n* + ^ /r2 ) % , 

D [W,W,W,W + c(i; 2 z + C* + ^ /r 2) % . 


C = T^T" 1 (9a) 

U - + E^+ 5j /r *> % ] 

V - c(n* + + n‘ /r >'" ] 

w - c( C | + c; /r > % ] 


(9b) 


with 
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^ ^ n, 

where 0 = k z u + k^v + k^^w and, for example, 

k z = k z / (k 2 + k 2 + k 2 /r *) % , etc. and k = ? for A, 
k = n for B, and k - 5 for C. Also <J> 2 = 0 .5(y-1) (u 2 +v 2 +w 2 ) , 
a = pc/2 and 3 = 1/a. 

The similarity transformations are applied to Eq. (8) 

and then a modified form is obtained by moving T^, T^, and 

T outside of the difference operators S_, 6 , and <5^, 

<= 5 n ? 

respectively. This results in the diagonal form 

T^ n 1 1 + hS^ n |N|l + h6 n A n n |P|l + hS^A^KT- 1 )^ 11 = R n (10a) 


with N = T" X T 

, N 1 = 

= T _1 T r , P 

= T _1 T , P _1 

rT 

1 

EH 

II 

5 1 

1 

n 5 

n c * 

q n ' 

where 





m 

1 

ra 2 

m 3 

-y m 4 

ym 4 

-m 2 

m l 

m 4 

ym 3 

-ym 3 

-1 





T T = -mo 
k z 3 

" m 4 

m l 

-yn^ 

ym2 

^ m 4 

-ym 3 


y 2 (1 + m^) 

y 2 (1 - m 1 ) 

-’P m 4 

ytn 3 

-yn^ 

y 2 (1 - m^) 

y 2 (1 + m^) 

m.. - k 1 

+ k 1 

+ k 1 , 

m_ = k 1 

k 1 , 

1 XX 

y y 

z z 

3 x z 

z X 

11 

K 

t-*! 

- k 1 
y x 

> 

m 4 = ' 

k_l , 
z y 

and y = 

1//2 





(10b) 


The source term and smoothing terms are dropped for 
simplicity . 
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The right-hand side, R , is exactly the same as in Eq. (8) 
and is formed in the same way. The implicit part of the 
algorithm now consists of 5x5 matrices, (i .e . , T^) and scalar 
tridiagonal operators, since, for example, |l + h6^A^| 
is made up of diagonal blocks. Therefore, the implicit part 
of the solution process consists of 3 matrix-vector multiplies 

A 

(i.e., T^ R) and 5 scalar tridiagonal inversions per direction 

(i.e., 1 1 + h6 | ~^T“^R) . This produces a substantial 

reduction in the total number of operators required for 

the implicit algorithm over the standard form. 

13 

Pulliam and Chaussee show that the diagonal algorithm 
retains the stability and convergence properties of the 
standard algorithm. Also, since the right-hand side is 
identical for the two algorithms, the steady-state solution 
is not affected by the diagonalization . Reductions in CPU 
time of up to 30% can be realized by using the diagonal form. 


4.3. Geometry and Mesh Generation Descriptions 

The generation of a computational mesh for calculating 
the flow through and around a Prop -Fan configuration consists 
of three steps. Step one is composed of generating an x-y-z 
system of points given the coordinates of the nacelle and 
blades. Step two consists of reading the x-y-z locations 
of step one and forming a system of parametric cubic patches. 
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Step three takes the patches of step two and based on the 
specified clustering and/or stretching of points, creates 
a mesh which is used by the computational code. 

The geometry for the blade, spinner, and nacelle are 
provided by the user. These characteristics, shown in 
Figures 3 and 4, consist of the blade twist in degrees, 

AB, based on the blade angle at 757* blade radius, the local 
blade chord ratio, b/D, and the leading edge alignment ratio, 
LEA/D, where D is the tip diameter. The thickness ratio 
distribution for each section is also needed in tabular 
form. The ultimate goal of this part of the mesh generation 
scheme was to smoothly transition from the twisting and 
tapering of the blade to a set of surfaces which have zero 
twist or taper above an r/D = 1 and with slope dcf>/dz = 0 
at the stagnation point of the nacelle and a zero slope at 
some specified distance downstream from the blade trailing 
edge. The two transition curves which are employed in this 
part of the analysis are parabolic and cubic curve fits. 

These curves allow for end points and/or slopes at the end 
points to be specified in the transitioning process. For 
the case of a parabola the following equations and conditions 
are used. 
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Figure 3. Prop-Fan 
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Figure 4. Prop-Fan Blade 
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General Equation: 

Y - aX 2 + bX + c 

Given : X. 

X 

initial 

final 

__ > 

, and 

initial 

Y final 


Final form of the coefficients : 


dY 

dX 

initial 


a 



- Y i> 

(X f - V* 


b 



(x t + X f ) 
CX i - X f ) 



< Y f - Y i> 
(X f - X ± ) 2 


( 11 ) 


c = Y f - x f y: 


X 


+ (X f - 2X i ) 


< Y f ~ Y j> 
(X- - X,.) 2 


L < x f - x i> v-f 

where subscripts i and f stand for initial and final conditions, 
respectively. 

If a cubic is used, the following equations and conditions 
are used. 

General Equation: X = aY 3 + bY 2 + cY + d 


Given : X. . ^ . . 

initial 


X final ^initial 


Y 

final 


X final 


initial 

Final form of the coefficients : 

2 (X f - X ± ) - (Y f - Y ± ) (Xf + X*) 


b = 


-< Y f ' Y i > 3 

- (X f - X i> - X i< Y f " Y i> 


0T f - Y.) 


- a * (Y- + 2Y. ) 


c = X^ - b • (2Y i ) - a • (3Y|) 


d = X f - a • (Y 3 ) - b • (Y 2 ) - 


00 


( 12 ) 
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where the subscripts i and f are the same as before. 

A user has to be careful when using cubics or parabolas. 

Once these surfaces shown in Figure 5 have been 
specified, it is an easy matter to extend the surfaces at 
constant values of any specified outer boundary. Usually 
this boundary is chosen to be far enough removed from the 
configuration such that there are no outer interactions 
effecting the computation. If a tunnel wall is to be specified, 
this would become the outer boundary. 

This whole procedure is automated with only outer boundary 
points and/or slopes specified in addition to the coordinates 
of the blade and nacelle surfaces. These parameters are 
presented in Figure 5. 

Step 2 is a self sufficient part of the grid generation 
program. An x-y-z set of points from step one is all that 
is necessary to create a series of parametric cubic patches 
which describes the entire mesh system in terms of cubic 
equations . 

In this section, an attempt will be made to briefly 
familiarize the user with terminology and definitions 
associated with the creation of a parametric cubic patch 
interpolation surface. This section is not intended to be 
a comprehensive narrative on parametric cubic interpolation. 

A far better approach to understanding the subject would 
begin with a thorough study of curves . Since this would be 
beyond the scope of this report, only three-dimensional 
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surfaces will be treated. The main purpose is to orient the 
user prior to presenting the user instructions in a following 
section. 

The basic element of the mesh generation code is the 
parametric bicubic patch. The patch represents a continuous 
mapping of a portion of the three-dimensional surface into 
a unit square by means of the equation 

4 4 

r(u,w) = ^ ^ ^ ^ A^. u 1 1 w^" 1 (0 <. u, w £ 1) (13) 

i=l j=l ^ 

This rather simple equation lends itself to fast, efficient 
computation. From a practical standpoint, the most important 
aspect of this equation is the determination of the coefficients 
A- . The coefficients A. are referred to as the algebraic 

j j 

form of the patch coefficients . Although this form permits 
efficient evaluation, it reveals very little information 
about the patch or how a user might determine the A. . Using 
matrix notation, Equation (13) can also be written in the form 

r(u,w) = UMBM T W T (14) 

o o o o 

where U and W are the vectors (u , u , u, 1) and (w , w , w, 1) , 
and the matrix, M, is given by 

2-211 
-3 3 -2 -1 

0 0 10 

10 0 0 
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The vector matrix S is referred to as the geometric form of 
the patch coefficients and consists of the following elements: 


*(0,0) 

*(0,1) 

!*(0,0) 
O w 


*(1,0) 

r(l.l) 


5£(i,i) 

3w 

f^CO.0) 

3u 

|i(o , 1) 

3u 


1 5 (o,D 

3u3w 

^(1,0) 

3u 


32? 0,0) 
3u3w 

8 * (1,D 
3u3w 


Most notable about the geometric form of the patch coefficients 
is that it contains information evaluated only at the patch 
comer points. In this form, all elements but the lower 
right-hand quadrant in Equation (16) become geometrically 
identifiable as comer point locations and tangent vectors. 

In fact, using the boundary curve numbering convention 
shown in Figure 6 the first two rows of B completely 
specify boundary curves 1 and 3 while the first two 
columns specify curves 2 and 4. 

Thus, 12 of the 16 elements of B may be determined from 
knowledge of the patch boundary curves. It should be mentioned 
that while Figure 6 shows the patch having four distinct 
boundary curves , degenerate cases often occur in practice 
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where one or more of the curves reduce to a point. The 
mixed derivatives control the behavior of the interpolation 
surface in the interior of the patch. When Coons 
originally formulated the patch concept, the mixed derivatives, 
or twist vectors as they are referred to in the literature, 
were identically zero. Thus , the entire patch depended only 
on the nature of the boundary curves. Although the simplicity 
of this concept is attractive, the attendant loss of capability 
is considered unacceptable. 
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One of the more subtle features of parametric interpolation 
concerns the parameters themselves. The unspecified relation- 
ship between the parameter and physical arc length creates an 
additional degree of freedom. To illustrate this point, 
consider a single curve segment to be represented by a 
parametric cubic of the form 


r(fi> = Y, 
i=l 



( 17 ) 


dr 

The coefficients A. may be determined by specifying r and — 

i du 

at the curve end points. At first glance, this appears to 

be equivalent to specifying the end-point locations and slopes. 

However, the end-point slopes or direction cosines are 

represented by unit vectors — ^ where s is the arc length 

ds 

along the curve . Thus , 


dr _ dr ds 

clu 3s" du 


( 18 ) 


An entire family of parametric cubic curves can be created 

which satisfies the same end-point conditions by varying the 

ds 

assumed value of — at either end-point. Specifying end-point 
locations and slopes for an ordinary cubic (e.g., y(x)) 
determines a single unique curve. 

The determination of the twist vectors at the patch 
comer points is much more difficult since they are not as 
geometrically identifiable as are the slopes or point locations. 
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For this reason, the methods used in the present program 
remove this burden from the user. 

Finally, this allows for full geometric derivative 
information available for any part of the mesh since the 
cubic functions are analytical in nature. The coefficients 
of the cubics are stored off line during the computational 
cycle to be available whenever needed in step 3. 

In the third and final step, the computational mesh 
used in the flow field code is created from the series of 
patches in step 2. 

Capsulizing, the overall procedure, steps one through 
three, has been automated (overlayed on the computer) to 
the extent that a user need only furnish the x-y-z points 
of: the blades, the location of the upstream, outer and 

downstream boundaries, and the spacing of the points in 
such critical regions as the leading and trailing edge of the 
blades, the tip region of the blade, and the region in close 
proximity to the nacelle surface. The final mesh is then 
created and stored on disc or tape for future use by the flow 
field code. It should be noted that the mesh should be 
displayed via computer generated plots. This allows for an 
evaluation of the grid system before it is actually used 
in conjunction with the flow field code. 
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4.4. Boundary and Initial Conditions 

The interior of the computational domain is solved 
implicitly by integrating the Euler equations. To complete 
the problem, boundary conditions must be specified on all 
surfaces of the computational volume. On the impermeable 
surfaces (blade, spinner, nacelle) , a surface tangency 
conditions is imposed. This consists of employing schemes, 
with as much accuracy as possible, such as Kentzer's method, 
a normal momentum approach as is done in the present code, 
or the method of characteristics. At the upstream and down- 
stream boundaries, anechoic like boundary conditions will 
be imposed. Free-stream values are used at radial far-stream 
boundaries with the exception being tunnel walls. These 
would require the tangency condition to be applied. In the 
region of the blade leading and trailing edges and at the 
blade tip a Kutta condition is used as in the usual airfoil 
analysis to allow for the flow to leave the surface of the 
blade smoothly. 

The following sections discuss ways of making the 
tangency condition at the blade and nacelle surfaces more 
exact, applying a three-dimensional Kutta-like condition 
and applying a subsonic outflow condition. 
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4.4.1. Tangency Condition 

Presently a straightforward approach is used at both the 
blade and nacelle surfaces to satisfy tangency. This consists 
in the case of the blade surface, of specifying the appropriate 
contravarient velocity to be zero (no flow through the surface) 
w = 0 and solving, along the £(x,y,z,t) = constant surface, 
the system of equations 


u 


^r n 4> n r^^ r 


5-5 t 

V 

= J" 1 

-< n zVV* )/r a z Wz )/r ■ <5 zW Z )/r 


V - n t 

w 


(n i ? r" n r' z ) ' (5 z ! r-Vz> 


w-c t 


(19) 


to obtain the cylindrical components of velocity u, v, and w. 

To solve this system, a complete set of contravariant velocities 
must be known and thus the u and v contravariant velocities 
are linearly extrapolated from known data on the interior. 

The pressure is obtained by zeroth-order extrapolation from 
the interior field and the density is calculated either by 
assuming constant entropy on the surface or by zeroth-order 
extrapolation. This completes the set of equations for the 
boundary conditions at the nacelle and blade surfaces. 

Unfortunately due to the extrapolation of the variables, 
the condition of total rothalpy equals a constant throughout 
the flow field is not satisfied. The total rothalpy is defined 
in a rotating system as 
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H = H - 2rww 
R t 


( 20 ) 


where =» -X— £ + l(u 2 + v 2 + w 2 ) is the definition of total 
u Y-l P Z 

enthalpy in a non-rotating system. It can be seen that 
rothalpy and enthalpy equal each other if the system is not 
rotating. In a uniform free stream, the total rothalpy is 
defined as 

h r = ^ + + ^ + w “ ) ( 21 > 

and it is held fixed at this value throughout the computational 
domain. In the following development, it is assumed the 
pressure was obtained by extrapolation, the density from the 
constant entropy relation 

1 

P = ( § )Y (22) 


where S = — 

Y 

Poo 

and the cylindrical velocities by solving the system of 
Eq . (19). There are three surfaces along which tangency 
has to be satisfied: 

I . Blade surface 

II. Nacelle surface 

III. Blade-nacelle intersection. 

At the blade surface, the contravariant velocity 

w = c t + us z + v? r + | ^ - 


0 


(23) 



39 


and the total rothalpy, H 

R 

u 2 + v 2 + w 2 - 2rtow = 20L, - 2-) (24) 

K y-1 p 


are known. Since there are two equations in the three 
unknowns u, v, and w, one of the velocity components, u, 
is accepted as being correct. The remaining components, 
v and w are found by solving Eqs . (23) and (24) , 
simultaneously. Doing this gives the v-velocity and 
w-velocity in terms of the accepted u-velocity and they 
are written as 


w = rco + -b ±yj b 2 - 4ac 

2a 


v = - 


il 


vlC, + (w - rw) 




] 


(25) 


where a = 


b = 


c = 


1 

2ur£ 


’* L 


2 t 2 




{2(H - -i- E) 

r K Y-1 P 


- u 2 + r 2 co 2 } 


■] 


By using Eq. (25) , the variables at the blade surface satisfy 
both surface tangency and the fact that total rothalpy equals 
a constant. 

At the nacelle surface, the contravariant velocity 


v = n t + u n z + vn r + ^ = ° 


(26) 
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and the total rothalpy, Eq. (21), are known. In this case, 
the w velocity is accepted and Eqs . (24) and (26) are solved 
simultaneously to yield 


and 


2rww - w 2 + 2 (H^ - 


_L_ R) 

Y-l P 


% 


u = 


1 + ^ 

nr 


(27) 


v = 



(28) 


Once again, the conditions of surface tangency and total 
rothalpy equals a constant are satisfied when Eqs. (27) and 
(28) are used in conjunction with the fact that w, p, and p 
are known. 

The final case occurs when the blade and nacelle 
intersect. At this junction both Eqs . (23) and (26) are 
satisfied due to a simultaneous condition of surface tangency 
on the blade and nacelle surfaces. Accepting the pressure, p, 
and the density, p, the cylindrical velocity components are 
calculated from 


v = 


( 2{H r ~ + r 2 ui 2 ) 


L 


z 

2 (n 2 


% 


r + n z } + r 2 (n r c z 


(29) 


Vz y 


and 



u 


(30) 
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These velocity components satisfy both surface tangency and 
the condition of total rothalpy equals constant at the blade 
nacelle intersection. 


4.4.2. Kutta Condition 

This is a hypothesis put forth by Kutta and Joukowsky 
to allow the inviscid flow past an airfoil to have a form 
approximating an ideal one and yet be physically possible. 

This hypothesis is necessary since no fluid can flow around 
a sharp trailing edge' — because any velocity around an arc 
of zero radius would mean an infinite acceleration. The 
flow no sooner starts than it finds itself sliding tangentially 
past the trailing edge which approximates the physical 
situation . 


C 

A - A 

(a) (b) 

Figure 7. Kutta Condition Geometry 
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In general, the unit normal to a 
Figure 7, is written 


n = 


+ ^r^r + r S 
K + ? r + { ^ }2)% 


? 


constant surface, 


(32) 


To implement the Kutta-Joukowsky condition, the unit normals 
of the face side and of the camber side are calculated and 
then these vectors are used to form a plane which is normal 
to the trailing edge of the blade and any given radial 
location. The resulting procedure is as follows. First 
form the unit normals 


n c - 


K + ! ^ } *> 


(33) 


camber 


and 



a y a 

‘’z^z 

ct‘ z + A 2 ) % 


face 


(34) 


Since two vectors from a plane, the scalar product of 
Eqs . (33) and (34) is performed to give 

1 

n c • n f = \zf c \ |^ f | cos 0 


(35) 


which defines the angle 0 between these vectors in the plane. 
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Figure 8. Vector Representation 

The desired result is finding a -vector which bisects the 
unit normals n c and n^. This is written 


1 


/V 

n c 

‘ cos 0 / 2 

(36) 


1 


A 

n f 

' = l^f| cos 6/2 

(37) 

A 

"5 

= 1 

(38) 


Equations (36), (37), and (38) form a set of three equations 
in three unknowns and are cast in the following form 


n % " 

ai 
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+ 

bi 
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Cl , 
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.a 
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/% 
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dl z 

+ 

ei r 
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fi 
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/\ 
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A 

ri 

z 

+ 

/N 

si 

r 

+ 

A 

ti 
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( 39 ) 
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where the coefficients d, e, f, r, s, and t are known from 
the metrics and a, b, and c are the unknown quantities. From 
Figure 8 the angle 0 is determined from 

n c -n^ = cos 0 = d*r + e-s + f-t (40) 

the system of Eqs. (39) and (40) are solved for 

a = cos0/2 - e-b - f-c 

b = A - c-B (41) 



it ^ d 

where A = ( — -) cos 9/2 

er - S • d 

B = ^‘ r ~ d 

e • r - S • d 

A c = d 2 (B 2 + 1) + f 2 + e 2 -B 2 - 2f *e -B 

B c = 2{(eB - f) cos 0/2 + e-f-A - e 2 'AB - A-B*d 2 } 

C c = d 2 (A 2 - 1) + e 2 A 2 + cos 2 0/2 - 2e-A-cos 0/2 

The velocity vector at the trailing edge is required to lie 

in the direction of the unit normal n, . The pressure is 
obtained by averaging the values from the camber and face 
side of the blade and the density is obtained from the 
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definition of total rothalpy after the velocity components 
have been adjusted such that they vectorially lie in the 
direction of the unit normal n, which bisects the unit 
normals of the camber and face side of the blade. 

This same procedure with slight modifications could be 
used at the tip of the blade to assure the flow leaves the 
tip smoothly rather than attempting an infinite acceleration 
around the tip. 


4.4.3. Subsonic Outflow 


The relationships enforcing boundary conditions at the 
subsonic outflow are determined from the governing partial 
differential equations using the theory of the method of 
characteristics (MOC) . The flow for this study is considered 
to be nonviscous and non -heat conducting. The equations cast 
in cylindrical coordinates for nonviscous, non -heat conducting, 
adiabatic, compressible flow are given below. 


Continuity 


l£ + + 1L 3P + P 3 v+v3£ + P3w,w 3£ = _p v 

3t P 3z 3z r 3r r x 3<f> r 3<j> r 
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z- Momentum 
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( 43 ) 
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d> -Momentum 
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(45) 


Energy 
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3z 

<h 
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(46) 


where c is the isentropic speed of sound for a perfect gas. 

Equations (42) - (46) are in cylindrical coordinates, 
thus in order to obtain these equations in a generalized 
coordinate frame of reference, the following transformation 
must be applied. 

t = t 

K = 5(t,r, <j>, z) 

n = n(t, r, <j>, z) (47) 

C = c (t , r, <f>, z) 

By applying transformation (47) , the partial derivatives in 
Eqs . (42) - (46) may be written generally as 
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( 48 ) 
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a 4» ax 3<j> 35 a<j> an a$ 


an , a 
— + ac 3<t> 


Combining Eq. (46) and Eq. (42), Eqs . (42) - (45) are 
written in generalized coordinates utilizing the derivative 
definitions of Eq. (48) . 


Continuity 
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where the contravariant velocities are defined as 


u = 


H + v ii + w 
+ + v 3r + r H 


v = 


w = 


fa + u|n + + 


,3n , W 3n_ 
3r r 3(j> 

3 1 3z 8r r 3 ^ 


( 53 ) 


The method of characteristics will now be formulated by using 
the set of Eqs . (49), (50), (51), and (52). A linear 
combination of the continuity Eq. (49), the 5 -momentum 
Eq. (50), and n -momentum Eq. (51) and the z, -momentum equation 
from (52) is obtained by multiplying them by scalars yp p£, 
1*3’ ^ 4 » respectively. The new vector coefficients of these 
equations now contain the scalar quantity pp The plane 
containing the new vector coefficients of the equation is 
called the characteristic plane. A normal to this plane 
is "A and is called a characteristic normal. The dot product 
of the new coefficient vector and the normal vector A must 
be zero. This procedure will be followed to determine in 
the characteristic compatibility equation in the 5-direction 
for the exit boundary condition. 

Rewriting Eqs. (49), (50), (51), and (52) 


w. 


P T + up 5 + pc 2 (u ? c z + v s £ r + — \ 


(54) 


w 


w 


■ v Pp - wPj - »' ! (V 2 + V r + 5T % + u C C z + v C«r + T 5, Cf) 


49 
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v T + 5v E + i p 5 5 r = r 3 = -w n - wv ? - I (p n n r + P ??r ) 
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(55) 

(56) 

(57) 


A characteristic like compatibility equation in physical space 
can be derived by making a linear combination of the above 
equations . 
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For a characteristic-like equation, the following relations 
must hold 

"i 5 ^2r +l, 3r + > 4l = >1J i 

y l pc2 ^z + y 2^ 

Vl pc2 C r + h 3 u 


= Xy 2 

(59) 
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or rewriting Eq. (59) 
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A nontrivial solution exists if the determinant of the 
coefficient matrix vanishes. The expansion of the coefficient 
matrix determinant gives the characteristic equation 


(u-A ) 2 


(u-A) 2 



q + q + 



= 0 


The roots of this characteristic equation are 
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The compatibility equation in the ^-direction is written as 


y i (p T + Ap ? } + p 2 {u x + X V + y 3 {v T + X V 

4 

+ y, {w + Aw } - E y. R. =0 (64) 

4 T 5 x i 

where A and y^ are defined in Eq. 


(63) . 
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For the outflow (exit) boundary condition, which is a 
finite distance downstream of the blade trailing edge, a 
constant pressure is specified. This specification is c 
consistent with the physical aspects of the problem and the 
theory of characteristics. The u-c characteristic brings 
information into the central volume from outside when u < c. 

Thus one dependent variable such as the back pressure must 
be specified. The remaining dependent variables are determined 
by using the up running characteristic (+) in the compatibility 
relationship, Eq. (64), to solve for the axial velocity, u. 

The term is evaluated using data from the previous time 
step. The momentum equations are used to solve for the 
remaining velocity components, v and w, and using the energy 
equation to solve for the density. 

The uprunning characteristic form of the compatibility 
equation, Eq. (64), with the appropriate A and y^'s is written as 
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(65) 
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With p specified as a constant in time, p T = 0, Eq. (65) 
is solved for u T 



( 66 ) 


For the grid being considered in this report, it is known 
that = 0 at the downstream boundary, causing a 

simplification of Eq. ( 66 ). If this was not the case, 

Eq. ( 66 ) would have to be solved simultaneously with the 
v- and w-momentum equation. 

Writing equations with the simplification, £ r = = 0 

the equation for the u-velocity is obtained. 

u T - - (u + CC Z )U E + *2 + ^ ( 9 / + « z P E > (67) 
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The v and w velocities are obtained from the momentum 
equations (51) and (52), 

v t = - (uv ? + w + vv ? ) - i (p ? ? r + p n n r + p ? 0 (68) 

and 


W T “ - (uw 5 + ™n + i + P n % + P;S >/r (69) 

The density is then determined from the following form of 
the energy equation 




(70a) 


where ± s the substantial derivative. Expanding Eq. (70a) 
and making the appropriate substitutions a characteristic- like 
equation for density is obtained, 

<P T + up ? + vp n + wp^) 

P T = — - - vp n - V/P c (70b) 

Since pressure, p, is specified, Eqs . (67), (68), (69), 
and (70b) are used to solve for the remaining independent 
variables u, v, w, and p at the subsonic exit plane. 


4.4.4. Initial Conditions 

There are two possibilities to be considered: 

1) Startup from a specified free stream, 

2) Startup from a previously calculated solution. 
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When starting from a specified free stream, the free 
stream variables are specified at every nodal point at time 
T = 0. As marching proceeds in time from T = 0, the tangency 
condition is imposed according to the cubic formula 

S T = 10 - 15T + 6T 2 + T 3 (71a) 

as applied to the velocity components 

u T-f.i - u co(l • 

v T+i = v T S T + v m ( 1.-S t ) (71b) 

w T+i = W T S T + 

with the following restriction imposed, T = 1.0 then S T - 1.0 
where 

T = (Iteration step) /(number of steps to slow (71c) 

start tangency B.C.) 

The remaining variables are accepted as predicted by the 
numerical algorithm in conjunction with the total rotholpy 
and entrapy requirement. 

This procedure is used over a finite number of time 
steps to allow for a smooth adjustment from a constant 
field to one with a large perturbation caused by the body 
imposed upon it. 
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When starting from a previously calculated solution 
using new specified flow conditions, the procedure is the 
same as before with the exception being the u^, v^, and w^ 
velocity components being replaced by the initial velocity 
components from the previous solution . 

If a solution is being continued in time with the same 
flow conditions, it is just a simple matter of reading all 
of the necessary information at the previous time level 
from the stored calculated solution. 

4.5. Output Data Manipulation 

Possibly the most important part of any design computer 
code is its output. To enable the designer to best evaluate 
a given configuration, the detailed surface pressures and the 
overall forces and moments on the blades and nacelle are 
necessary. The blade section aerodynamic coefficients and 
the overall performance coefficients of the blading are also 
of great importance in making the right choice of the most 
efficient airfoil section. In addition, computer generated 
plots of the various aspects of the flow field surrounding 
the Prop-Fan configuration would be beneficial in the design 
process . 
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Figure 9 . Surface Definition 
Since force and moment coefficients are required for both 
the blades and the nacelle, a generalized derivation is 
included in this section. The surface definition in generalized 
coordinates for the calculation of forces and moments is shown 
in the Figure above. The incremental force acting on the area 
segment of the surface, dA can be written 

dF = -Ap n dA (72) 

where the unit normal to the area segment is defined as 
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AAA 



The incremental area in generalized coordinates is 


(73) 



[ s 


i + s 2 + 
x y 


] 


d^dn 


(74) 


where J is the Jacobian. 

Thus the incremental force, dF, can be written in 
generalized coordinates as 


dF = - ^ (S i + S 
J x 


7 * 


S k) 
z 


d£dn 


(75) 


In order to determine the various components of the force, 
the scalar product incremental force dF with the unit vector 

AAA 

components, i, j, k, has to be formed. The usual convention 
for assigning forces is shown in the following Figure. 


z,N 



Figure 10. Force Assignment Conventions 
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The various components, N, S, A, of the total force F 
are defined as the normal force, side force, and axial force, 
respectively. To arrive at a mathematical definition of the 
various components, the scalar product of the generalized 
incremental force dF with the unit vectors i, j, k is formed 
to give 

dN = dF-k = - £2. s z d£dn 

dA = dF'i = - d^dn (76) 

J x 

dS = dF* j = - AE d£dn 
J j y 

and integrating these incremental forces, the total force in 
each direction is obtained as follows 


Normal Force 


N = 

Axial Force 

A = 

Side Force 

S = 


Iff S Z 

n K 

Iff S x d ^ dn 

n 5 


Iff 



d£dn 


(77) 


(78) 


(79) 


The incremental moments dM are obtained by forming the 

4 1 

vector product of the radius r with the incremental force 
dF to give 
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dM = r x dF (80) 

where 

-+ S' Ss * 

r = xi + yj + rk 


Thus, the following equation is obtained for the total 
incremental moment, dM, 

dM = - (yS z - zSy)i + (zS x - xS z ) j + (xS y - yS x )k] d£d n 

(81) 

and as in the forces, the various components of the total 
moment can be written as 

dM^ = dM’i 

dMp = dM-j (82) 

dMy = dM'k 

where the sign convention for these moments is defined in 
Figure 4. Integrating these incremental moments yields 
the following total moments measured with respect to the 
nose of the body 

Rolling Moment 

Mr - - fff- (yS 2 - zs y )d 5 dn (83) 

K n 

Pitching Moment 

«p - - fff- (zS x - * s z> d «*> 

S n 


( 84 ) 
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Yawing Moment 


My “ - Jf f- (.-xSy - yS x )d5dn (85) 

5 n 

Typically, the center of pressure may also be of interest 
and it is defined as 

v = PITCHING MOMENT ^p (Q 

C P NORMAL FORCE ‘ N Q ' 


Also, the forces and moments are converted to nondimensional 
coefficients as follows 


C 


F 


Force 

(Dynamic Pressure) (Reference Area) 


°M ~ 


Moment 

(Dynamic Pressure) (Reference Area) (Reference Length) 


(87) 


where 


1 2 , 

q =2 Pco v «> = dynamic pressure 
A = iTr^ = reference area 
= reference length 


The various forces and moments are then written in 


coefficient form as 
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<lo A b 

_N 

%> A b 


_S 

‘leA 


( 88 ) 


and 



E_ 

‘UA/b 

\ 

q ~Vb 

My 

q<,Ab £ b 




(89) 


This generalized derivation of the force and moment 
coefficients will be applied in the sections to give values 
for the blades and nacelle. 


4.5.1. Blade Properties 

The blade properties are broken up into three distinct 
categories, each giving valuable information to the designer. 
In the first category, the aerodynamic coefficients of the 
blading are calculated. These coefficients consist of the 
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axial force, normal force, and side force along with the 
pitching moment, rolling moment, and the yawing moment. 

The second category defines the performance coefficient, 
torque coefficient, and the power coefficient. The 
independent variables compose the three category. These 
are generally the pressure, density, velocity, and energy 
variations on the blade surface. This category would be 
considered as a standard output with categories one and 
two being more specifically aimed at the Prop-Fan problem. 
With this in mind, a description of both category one and 
category two follows . The geometric relationships are 
shown in Figure 11 below. 



Figure 11. Geometric Relationships for 

Conversion to Cylindrical Coordinates 
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The generalized cartesian coordinate development of 
forces and moments is now converted to cylindrical coordinates 
using the relations 

x = z 

y = R cos<j> (90) 

z' = R sin<}> 

and the fact that the cylindrical Jacobian is equal to R times 
the cartesian Jacobian, 

J cy x - R-J < 91 > 

Using the coordinate transformation Eq. (72) the 
cylindrical Jacobian Eq. (91) and the generalized forces 
and moments Eq. (76) and (82), the incremental forces and 
moments in cylindrical coordinates are 

dF AXIAL “ ' J^ d?dn 
dF N0RMAL = ' (cos <|>S R - 


dF SIDE " j P (sincj)5 R + ? ^d^dn 

cyl *'• 


^ROLL " d?dn 
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(z{cos<j>S R - - Rcos<K z >d£dn 


^AW 



(-z{sin<j)5 R + 


cos 4> 
R 



+ Rsin<j)£ z ) d£dn (92) 


The equation set (92) is integrated in space by using 
the trapezoidal rule to obtain the total force and moment 
coefficients contributed by the blades. 

Propeller characteristics are normally defined in terms 
of efficiency, n, thrust coefficient, C^,, torque coefficient, 
Cq, power coefficient, Cp, and advance ratio V/(nD) . The 
determination of these characteristics is presented below. 
The notation is given in the following Figure. 



Figure 12. Notation for Propeller Characteristics 
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In the Figure 

6 = geometric blade angle, based on chord, at station r 
cf> = effective pitch angle 

a = effective angle of attack = (90° - g) - $ 

V = free-stream velocity 

fir = linear velocity due to rotation at station r 

dL = elemental lift 

dN = elemental normal force 

dD = elemental drag 

dA = elemental axial force 


From an examination of the Figure, it will be noted that 

the elemental lift and drag forces are taken respectively 

perpendicular and parallel to the velocity vector v . Because 

R 

the velocity V R is the resultant of V and fir, it represents 
the net flow past the element. 

The elemental thrust is the force produced by the blade 
element along the line of flight, or 

dT = dLcos<J) - dDsincf) ( 93 ) 

The elemental lift and drag can be written in terms of the 
normal force and axial force as 


and 


dL = dNcosa - dAs in a 


(94) 


dD = dNsina + dAcosa 


(95) 
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This form of the coefficients is chosen since the forces and 
moments in this report are derived to be consistent with a 
body axis system rather than a wind axis system. Therefore, 
the elemental thrust is 

dT = dNcos(a+<j>) - dAsin(a+<j>) . (96) 

The elemental torque is the force produced by the blade 
element resisting rotation multiplied by the radial distance 
r to the center of rotation, or 

dQ = r dLsincj) + r dDcos<j> (97) 

Using Eqs . (94) and (95) in Eq. (97), the elemental 
torque is written as 

dQ = r dNsin(a+<j>) + r dAcos (a+cj>) (98) 

where the elemental normal force and axial force are defined 
in Eq. (92) . Also it may be seen from Figure 12 that 

V = = ^5- (99) 

R cos<p cos? 

and 

<p = tan -1 t~ = tan" 1 (100) 

2Tir 


where 


D 


2R — blade diameter 
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By definition. 


elemental thrust coefficient = dGj = 


dT 

pn 2 D 4 


( 101 ) 


elemental thrque coefficient = dC n = — ^ — (102) 

Q pn 2 I> 5 

where p is density and n is the frequency. 

The section thrust and torque coefficients are obtained 
by dividing both sides of Eqs . (96) and (98) by an incremental 
radius, dr, 


— = cos(a+(J>) - -j— sin(a+<j>) 
dr dr dr 


(103) 


and 


r M sin(a+<|>) + r sin(a+<f>) (104) 

dr dr dr 

Values of the overall thrust and torque coefficients Gj, 
and Cq are determined by a trapezoidal integration with 
respect to r of Eqs. (101) and (102) with the help of 
Eqs. (103) and (104). 

The efficiency of the propeller, n, is defined as 

= Power Output _ (Thrust) (Velocity) _ TV 
Power Input (27m) (Torque) 2?mQ 

FromEq. (101) and (102) 

T = C T pn 2 D 4 
Q = C Q pn 2 D 5 


( 105 ) 
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where 


Cj = thrust coefficient = 


Cq = torque coefficient = 


D/2 



o 


dCrr, 

— - dr 
dr 



d(k 

— 9- dr 
dr 


such that 



2uCq 


where J = — 
nD 

If we define a power coefficient, 



then 


C = power coefficient = &P^ r - = — 2 — 
P pn 3 D 5 pn 3 D 5 




J 


where from Eqs . (106) and (108) 


(106) 


(107) 


(108) 


C p - 2 ttC q (109) 

dGj, dC 

A typical plot of and — P- against r/D is shown 

dr dr 

in Figure 13 and is compared with the theory of Rohrback. These 
curves are integrated in r to give total thrust and torque 
coefficients. Figure 14 presented the calculated performance 
coefficients of the blades. 



Thrust coefficient, d(C T )/d(r) 


M =o = °- 8 37.5 SHP/D 2 J =3.06 

35,000 ft ISA 800 ft/sec tip speed Cp = 1.7 



0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 


r = r/D 


Figure 13. Thrust and Power Coefficient Loading as a Function of Blade Radius 
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properties of the given blading 


* * 

* GEOPETPIC PROPERTIES t 

* * 


NO, CP BLADES f 8 
BLADE DIAHETER t t.COOO 


*••#****«************«** 
* * 

* FL'ChF IELO PROPERTIES A 

• * 


HACh :MjHBER ,8900000 

PRESSURE •* 1,0000000 

DENSITY * 1,0000000 

VELOCITY * ,?a<»5728 

ADVANCE RATIO * 3,0600000 

THRUST COEFFICIENT * ",7906009 

TORQUE COEFFICIENT ? ,*537881 

PCKER' COEFFICIENT s J.R795S30 

BLADE EFFICIENCY b ,7023082 


Figure 14. Computer Output of Performance Coefficients 
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4.5.2. Nacelle Properties 

The independent variables, pressure, density, velocity 
components, energy, and the force and moment coefficients are 
the properties of interest in this section. Again, as in the 
blade properties, the independent variables are a standard 
output and will not be discussed here. The forces and moments 
are presented in terms of the cylindrical coordinates, Eq. (90) . 
The incremental forces and moments in cylindrical coordinates 
for the nacelle are 



To obtain the total force and moments that are contributed 
by the nacelle equation Eqs . (110) are integrated using the 
trapezoidal rule. 
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5. MESH GENERATION PROGRAM 

Before beginning the user's instructions, a few significant 
characteristics of the program should be described. In its 
present form, the program may be executed on a stand-alone 
basis. It has four OVERLAY links, the blade, nacelle, and 
grid boundary geometry link, geometry generation link, data 
card generation link, and the grid generation link with the 
last link requiring the most storage. At present, input for 
link I and for link II are stored on disc and tape and only 
changed when deemed necessary, i.e. completely new geometry 
definitions. In the following sections, a more detailed 
description of different facets of the program complex will 
be presented. 

5.1. Program Organization and Interaction 

This program is OVERLAYed with the four links being 
governed by a main driving program. This program calls each 
OVERLAY as needed and these occur in successive order of 
OVERLAY 01 through OVERLAY 04. The organization is such that 
the first link reads the blade and nacelle geometry data 
along with the definition of the outer boundary of the grid 
system. Thus in effect the complete outer and inner boundaries 
are defined here. These data are assumed to be fairly constant 
for a given problem and are stored on disc or tape (TAPE5) , 
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to avoid the problem of reading in large amounts of data from 
cards. The first link then creates a series of patches which 
link the outer boundary and the inner boundary. This information 
is stored on TAPE10 and is used in the second link. 

From TAPE 4, the governing parameters of link two are 
read which are used to create patch coefficients from the 
data which was furnished by link one . These parameters are 
assumed to be the same for most geometries which this program 
will be associated with. Thus, it has been stored permanetly 
to keep from reading an unnecessary number of cards. This 
link in turn stores the data for use in link three. 

Link three creates data cards on disc (TAPE3) which are 
read by link four. These data cards determine the number of 
points in the grid, the clustering of the grid lines, and 
other pertinent information in a grid generation process . 

Link four reads the data cards from link three which 
determine the size and location of grid points and reads the 
parametric cubic coefficient data from link two which describe 
the general three-dimensional domain and finally combines the 
two to determine the finite difference mesh which is used by 
the 3-D flow field code. The final system is written to 
TAPE 11. 
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OVERLAY 
0.0 


Executive Link 


OVERLAY 

01 


Inner and 
outer 
boundary 
link 


OVERLAY 

02 


Parametric 
cubic coef- 
ficient link 


OVERLAY 

03 


No. of points 
and clustering 
link 


OVERLAY 

04 

Mesh 

creation 

link 


Read TAPE5 



Read TAPE 4 


Governing P . C , 
parameters 
which remain 
unchanged 


Read initial 
clustering 
data from 
cards and no. 
of points 


Read the P . C. 
coefficients 
from mass 
storage 
as needed 


Creates patches 
between inner 
and outer 
boundaries 


Stores these 
patches as x,y,s 
points of each 
patch on TAPE 10 


Read TAPE 10 
Basic x,y,z 
points used 
to create P.C. 
coefficients 


Store these 
in mass 
storage 


Calculates new 
location of 
x,y,z points on 
body surfaces 
based on 
clustering 
information 


Store these new 
x, y , z, locations 
on TAPE3 for 
link 4 


Read TAPE 3 


Actual 
clustered 
x,y, z body 
coordinates 
used to 
create the 
grid 


Calculates 

finite 

difference 

grid 


S tores grid 
on TAPE 11 
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5.3. Program Input 

The program input is described here. A dictionary of 
the input variables is provided followed by a description 
of the different input data formats . The user should only 
have to be familiar with the input variables of TAPE2 and 
TAPE5 with the input variables of TAPE4 furnished once. 

5.3.1. Dictionary of Input Variables 

For a general description of the variables read from 
TAPE4, the user is referred to Reference 15 which was produced 
for NASA/Ames Research Center, Moffett Field, CA with the 
cognizant person being Reese Sorenson. The following list 
describes all input variables to the program using TAPE2 or 
TAPE5. Those for TAPE5 will be discussed first. 


VARIABLE 


DESCRIPTION 


AOD the axial location of the nacelle stagnation 

point with respect to the coordinate system 
origin, < 0; normalized by the blade diameter. 

BB(I) chord of the blade at each radial station. 

BETA34 twist of the blade at 75% span. 

the axial location of the boundary between 
patches 9 and 10 at some station downstream 
of the blades, > 0; normalized by the blade 
diameter. 


BOD 


I 


CCP 

COD 

D 

DBET(I) 

DDDP 

DOD 

EOD 

LEA(I) 

NB 

NBU 

NN 

NP 

NT 

RB(I) 

RN(N) 

XOB (N , K) 

YCOB(N,K) 

YFOB (N,K) 
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the radial location of the outer boundary; 
normalized by the blade diameter. 

the radial location where Patch 3 (Figure 2) 
becomes normal to the axial coordinate; 
normalized by the blade diameter. 

physical blade diameter. 

twist of the blade at each radial station. 

the axial location of the far upstream boundary 
< 0; normalized by the blade diameter. 

the axial location of the immediate upstream 
boundary < 0; normalized by the blade diameter. 

the axial location of the far downstream 
boundary; normalized by the blade diameter. 

blade axial leading edge location with respect 
to the coordinate orgin. 

number of planes in the radial direction defining 
the blade geometry. 

additional number of points between the blade 
tip (NB) and the outer boundary. 

number of points in the axial direction defining 
the nacelle geometry. 

number of points in the chordwise direction 
defining the blade geometry. 

total number of points in radial direction; NT = NB+NBU 

radial location of the BB(I) and DBET(I) . 

cylindrical radius of the nacelle at each axial 
station . 

axial location of the YCOB(N.K) and YF0B(N,K) 
coordinates; normalized by the chord BB(N) . 

radial coordinate of the blade camber side at 
XOB(N,K); normalized by the chord BB(N) . 

radial coordinate of the blade face side at 
X0B(N,K) ; normalized by the chord BB(N) . 

axial stations for which RN(N) is defined. 


ZN (N) 
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TAPE 2 input 
VARIABLES 

CFSL 

DETNAC 

DETTIP 

DXILE 

DXISL 

DXITE 

DZERAT 

JCONE 

JLE 

JMAX 

JSL 

JTE 

KMAX 

KTIP 

LMAX 


variables are discussed next. 

DESCRIPTION 


fraction of the chord at which "midchord" 
clustering occurs . 

actual physical spacing desired radially at 
the nacelle surface. 

actual physical spacing desired radially at 
the blade tip. 

actual physical spacing desired axially at 
the blade leading edge. 

actual physical spacing desired axially at 
the blade "mid-chord". 

actual physical spacing desired axially at 
the blade trailing edge. 

actual physical spacing desired meridionally 
at blade surfaces. 

number of cones (including z-axis) at nose 
of nacelle. 

integer location of blade leading edge; 
usually 11. 

number of mesh points in ^-direction; £ 45. 

= 0, no clustering at mid-chord; 

= 1, clustering at mid- chord. 

integer location of blade trailing edge; 
usually 27. 

number of mesh points in n- direction; <.21 
integer location of blade tip; usually 11 
number of mesh points in c-direction; <. 11. 
number of blades. 


NBLADE 
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PRT10 logical variable for OVERLAY 1 print option; 

.TRUE, print; .FALSE, no print. 

PRT20 logical variable for OVERLAY 2 print option; 

same as PRT10 . 

PRT30 logical variable for OVERLAY 3 print option, 

same as PRT10. 

PRT40 logical variable for OVERLAY 4 print option; 

same as PRT10 . 

5.3.2. Input Data Format 

The data could all be input from cards but due to the 
bulk of this data it was deemed best to store that data which 
does not change much on TAPE4 and TAPE 5 . The data on TAPE4 
governs the parametric cubic patch program and would only- 
change if a user had extensive experience with that part of 
the program. The data on TAPES is the blade and nacelle 
geometry definitions along with the axial and radial boundary 
information. This data usually remains constant for a given 
nacelle blade configuration. Consequently, only the clustering 
information is read from cards (TAPE2) . The format for the 
cards to be stored on TAPE5 are: 



CARD NO . 1 


Column no. 

1 


D 

CARD NO. 

(1615) 




Column no. 

1+5 



CARD NO. 3 

(2F10.5) 

There are NN of this 

type card. 

Column no. 

1 

11 




ZN(N) 

RN (N) 



CARD NO. 4 

(1615) 




Column no . 

1+5 

6 - 10 

11 - 15 

16 - 20 


NB 

NP 

NBU 

NT 

CARD NO. 5 

8F10.5 

There are 

NB of this 

type card. 

Column no. 

1 

11 

21 

~3T 51 


B32 

DBET(I) 

LEA(I) 

BB(I) 


CARD NO. 6 

(10X, 2 (F10. 4, 

2F10.6)) 

There 

are (NP/2)X(K) 

of this type card. 



Column no. 

1 11 

~2i 

31 


^1 

51 

61 


XOB (N ,K) 

YCOB(N.K) 

YFOB (N ,K) 

XOB (N+l ,K) YCOB(N+l,K) 

YFOB (N+l, K) 

CARD NO. 7 

8F10.5 








Column no. 

1 11 

21 

31 

41 

51 

61 71 




BET A3 4 AOD 

BOD 

COD 

DOD 

EOD 

CCP DDDP 




This concludes the card types which are to be stored on TAPE5 . 









Next the cards read from TAPE2 are presented 


CARD NO. 1 (4L5) 


Column no. 

1 + 5 

6 - 10 

11 - 15 

"""16 -" 

20 




PRT10 

PRT20 

PR130 

PRT 40 




CARD NO. 2 

(1615) 

• 






Column no. 

1 + 5 

6-10 11-15 

16-20 21-25 

26-30 

31-35 

36-40 

41-45 80 


JCONE 

JLE JSL 

JTE JMAX 

KTIP 

KMAX 

LMAX 

NBLADE 


CARD NO. 3 


1 

11 

21 

31 

41 

51 

61 

80 

DXILE 

DXISL 

DXITE 

CFSL 

DETNAC 

DETTIP 

DZERAT 



Column no. 
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5.4. Program Output 

Due to the nature and complexity of grid generating 
programs there are a large number of write statements in all 
four OVERLAYS. There are far too many to discuss in this 
manual. However, it is suggested that the user trigger the 
print options in all four links to get a feeling and under- 
standing of what output data is available. A shortened 
version of a print from OVERLAY 4 can be found in 
Appendix A. 3. 

5.5. Operating Procedures and Trouble Shooting Hints 

It would enhance the user's understanding of the validity 
of the grid system by displaying the grid using the local 
system plot packages. A detailed picture of the grid system 
can either ensure the grid or allow the user to quickly 
trace down possible problems. 

Typical things to check in case of problems are: 

• the coordinates of the blade and nacelle. 

• location of the boundaries with respect 
to the origin of the coordinate system. 

• the values of the clustering parameters 
read from cards . 

Any or all these could cause problems in the grid generation 
phase which is probably the most important part of the 
calculation procedure. 
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6. THREE-DIMENSIONAL PROP -FAN FLOW FIELD PROGRAM 

This program numerically calculates the flow field solution 
surrounding the Prop -Fan configuration much as a wind tunnel 
test calculates the experimental results. The benefit here 
is in the fact that a lot more detailed data is available 
for analysis from the present numerical results. These results 
consist of the density, velocity components, and energy which 
can be manipulated to determine any other necessary flow 
variable of interest. The following sections will present 
a more detailed description of the program. 

6.1. Program Organization and Interaction 

The organization is very logical and done in a straight 
forward matter. Usually, the mesh from the mesh generation 
program and a previously stored solution are attached in the 
starting process. The set of equations are then iterated 
in time to a steady-state thus yielding a wealth of information 
about the flow field. This solution is stored on disc (or tape) 
and can later be attached by either the flow field program 
or by the data reduction program. This allows the user at 
any given time to view the performance of the blade-nacelle 
configuration by using the flow field variables which were 
calculated from the 3-D flow field code. 


8 









I 
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6.3. Program Input 

The input required by the 3-D flow field program is 
described in this section. A dictionary of the input variables 
is provided followed by a description of the input data format. 

6.3.1. Dictionary of Input Variables 

The following list includes all of the variables which 
are input to the program. 

VARIABLE DESCRIPTION 

ADVR negative of the advance ratio of the blading, J, 

where positive is considered to being clockwise 
direction . 

BETA34 geometric angle, Beta, at the 3/4 span of the 

blade B 3 / i» . 

BLDIAM physical blade diameter, D^. 

CN Courant Number, used to determine marching 

s tep size. 

FSMACH free stream Mach number, M*, 

GAMMA ratio of specific heats, y 

IBC = 0, applies blade boundary condition; 

= 1, nacelle only. 

IGRID = 0, rescale by new Jacobian; 

= 1, old grid system. 

IMETHD has no effect in current version. 

IMETOT = 0, do nothing; 

= 1, print out metrics and exit. 

IPTCH = 0 \ grid is read from Tape 11 

=1/ by user in GRID. 
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IREAD 

IWRIT 

JLE 

JMAX 

JTE 

KMAX 

KTIP 

LMAX 

METH 

NBLADE 

NC 

NCHNGE 

NMAX 

PER 

PINF 

Q 

RINF 

SMU 

SMUIMP 

TAU 


= 0, start from scratch; 

= 1, read starting solution from Tape 1. 

= 0, nothing; 

= 1, store solution on Tape 2. 

integer location of blade leading edge; 
usually 11 . 

number of mesh points in ^-direction; <_ 45 . 

integer location of blade trailing edge; 
usually 27 . 

number of mesh points in p -direction; < 21. 

integer location of blade tip; usually 11. 

number of mesh points in ^-direction; £ 11. 

not used. 

number of blades. 

parameter read from Tape 1. 

not used. 

number of iterations to be performed. 

angle between blades in degrees. 

free stream pressure normalized by yPo,; = 1 / y 

3-D storage array of independent conservative 
variables read from Tape 1. 

free stream density normalized by p^; = 1. 

explicit smoothing coefficient; 0(At) . 

implicit smoothing coefficient; 0(2*SMU) . 

current value of time . 
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6.3.2. Input Data Format 

The data is either input from cards (Tape 5) or from 
disc (Tape 1 and/or Tape 11) . In Appendix A. 2 an example of 
a typical input is provided. In general, the input cards 
are as follows : 

CARD NO. 1 
Column no. 


CARD NO. 2 (615, 3F10.5) 


Column no. 

1 + 5 

6-10 

11 - 15 

16 - 20 

21 - 25 


IRE AD 

IWRIT 

IPTCH 

NMAX 

METH 


Column no . 

26 - 30 

313 

41 

51 80 


IMETHD 

CN 

SMU 

SMUIMP 


CARD NO. 3 
Column no. 


If IPTCH <. 0, the following card is read. 
CARD NO. 4 (715, 3F10.5) 


Column no. 

1+5 

6-10 

11 - 15 

16 - 20 

21 - 25 


JMAX 

KMAX 

LMAX 

JLE 

JTE 


Column no. 

26 - 30 

31 - 35 

36 

46 

TS 80 


KTIP 

NBLADE 

PER 

BETA34 

BLDIAM 


On the above cards, remember that all integers are 
right adjusted. 


(615, 3F10.5) 


1 + 5 

6-10 

11 - 15 

16 - 20 

80 

IBC 

IMETOT 

NCHNGE 

IGRID 



(5F10.5) 


1 

11 

21 

31 

41 

80 

FSMACH 

GAMMA 

PINF 

RINF 

ADVR 



T 
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The possibility of reading the input data from discs 
is a real and distinct option within the 3-D flow field code. 

If IPTCH > 0, the following data is read from TAPE11. 
From INPUT 

TAPE 11; JMAX, KMAX, LMAX, JLE , JTE , KTIP , NBLADE , 

PER, BET A3 4, BLDIAM 

and from GRID 

TAPE 11; X(KL, J) , Y(KL,J), Z(KL,J) 
are read. This data is used in employing and determining 
the computational grid. 

If the solution is continued from a previous solution 
the starting flow field is read from TAPE1 in INPUT as, 

TAPE1; Q, TAU, NC . 

6.4. Program Output 

In the following section, a brief description of the 
3-D flow field program output is presented. Typically the 
output can be broken up into various sections . The first 
section consists of printing the input data. This defines 
all of the options and the pertinant variables which are 
to be used during execution. Also, identifiers from the 
GRID subroutine and from the EIGEN subroutine acknowledge 
the reading of the grid and the determining of the stepsize 
respectively . 
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During execution (or marching) , there is a different 
set of data printed. This can be considered the second 
section. It consists of the L2 and maximum residual with 
its location, the step number and the normalized density 
at the body one point away from the stagnation point, the 
forces and moments on the blade and the nacelle, and finally 
printer plots of the various independent variables. The 
forces and moments and the printer plots occur every IFREQ 
iterations where IFREQ is set in the main program. 

The third section occurs at the end of the iteration 
and consists of a quite detailed prirt out of variables 
which are used to determine the boundary conditions on the 
blades and nacelle. This allows a user to observe these 
variables and to get a better understanding of the flow 
condition on the body surfaces . This is also useful in the 
process of debugging. 

Finally, the flow field is printed in the fourth section 
in two parts. The first part prints the entire flow field 
at the periodic planes and it prints the body data only at 
the interior planes . The output consists of the normalized 
independent variables p, p, u, v, w, e, M, the total enthalpy 
ratio Hp/H^ and the geometric descriptions Z, R, and <J>. 

The second part prints these same variables on the camber 
and face side blade surface. 
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6.5. Operating Procedures and Trouble Shooting Hints 

The only comment that should be made here on operating 
procedures is that if problems occur when starting the 
solution from a free-stream guess, the nacelle only should 
be calculated initially then allow the blades to enter into 
the solution at a later time. This makes for an easier 
transition of the flow field. 

Usually problems associated with the code stopping are 
well documented by an appropriate output. In the past, the 
major cause of a stoppage has been too high a Courant Number 
or not enough grid resolution at the tip of the blade. The 
Courant Number problem is easy to correct but the grid 
resolution is not unless one has an increased storage capab- 
ility. If increased storage is available, then adding more 
radial points at the surfaces would also enhance the chances 
for a better solution. 
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7. DATA REDUCTION PROGRAM 

This program is designed to take the flow field data 
produced by the Three-Dimensional Prop-Fan Flow Field Program 
and generate printer plots, store data for computer generated 
plots, and to calculate the various propeller performance 
coefficients. The following sections will describe these 
tasks in more detail. 

7.1. Program Organization and Interaction 

As is the case in most computer programs, the input data 
is read from cards, disc, and/or tapes before any calculations 
begin. This input data is that which has been furnished by 
both the flow field program (Section 6) and the mesh generation 
program (Section 5) . Once this interaction is complete, the 
data reduction program proceeds to manipulate the input data 
for the user. This results in: 

• printer plots of the C distribution on the face and 
camber side of the ^ blades at various specified 
radial locations; 

• printer plots of the axial Cp distribution on the 
nacelle surface at various ^meridional angles; 

• section properties on the face and camber side of 
the blades; 

• geometric parameters for each section; 

• integrated section properties to provide section 
left, drag, thrusts, and torques; 

• integrated section thrusts and torques to calculate 
the power coefficient and efficiency; 


T 


T 
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• the output of the above calculated pertinent quantities; 

• the three-dimensional data base being converted into 
various two-dimensional data bases which can be us ed 
by a plot routine to display the distribution of the 
various properties . 


7.2. Flow Chart 

This section provides a brief flow chart describing .the 
overall logical flow of the program. 
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7.3. Program Input 

In this section, the input required by the program is 
described. First a dictionary of the input variables is 
provided followed by a description of the input data format 
including an example. 

7.3.1. Dictionary of Input Variables. 


The variables which are input to this program are described 
in the following list . 


VARIABLE 


DESCRIPTION 


ADVR 

BETA34 

BLDIAM 

FSMACH 

GAMMA 

IPLOT 

JLE 


negative of the advance ratio of the 
blading, J. See Section 6.3.1. 

geometric angle, Beta, at the 3/4 span of 
the blade, B 3 / 4 . See Section 6.3.1. 

physical blade diameter, D^. 

free stream Mach number, M^,. 

ratio of specific heats, y 

control option for storing RSEC(N) 
on tape 9; equal to 0 or 1 . 

integer location of blade leading edge. 


JMAX number of mesh points in £ direction; <. 45 

JTE integer location of blade trailing edge. 


KLMAX number of KMAX*LMAX points; <. 240. 

KMAX number of mesh points in q direction; <_ 21. 
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KTIP 
LMAX 
N BLADE 
NC 

NMAX 

PERD 

PINF 

Q 

RINF 

RSEC(N) 

TAU 

X(KL,J) 
Y(KL, J) 
Z (KL , J ) 


integer location of blade tip. 

number of mesh points in £ direction; <, 11. 

number of blades . 

parameter read from tape 1; redefined within 
program. 

number of radial blade sections. 

angle between blades in degrees. 

free stream pressure normalized by YP^; = 1/Y . 

three dimensional storage array of independent 
conservative variables read from tape 1. 

free stream density normalized by p ro ; = 1. 

radial locations on blade for which flow properties 

will be printed; < 20. 

time the blade has traveled; not used. 

vector array of length KLMAX by JMAX of user 
supplied X-mesh points; read from tape 11. 

vector array of length KLMAX by JMAX of user 
supplied Y-mesh points; read from tape 11. 

vector array of length KLMAX by JMAX of user 
supplied Z-mesh points; read from tape 11. 


7.3.2. Input Data Format 


The data are input from cards (TAPE5) and from disc 
(TAPE1 and TAPE11) . An example of the input is provided 
in Appendix A. 3. The input data cards would read: 


CARD NO. 1 (8F10.6) 


Column no. 


T 

IT 

~n 

"31 

41 

80 

FSMACH 

GAMMA 

PINF 

RINF 

ADVR 
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CARD NO. 2 
Column no. 


CARD NO. 3 


Column no . 


(1615) 


1 - 5 

6-10 


80 


Max 

iPLOt 




(8F10.6) 

There are 

(NMAX-l)/8 + 1 

cards , 

i.e., NMAX = 10, 2 

cards . 


1 

H 

~ZI - - - 

zr ~71 go - 

RSEC (1) 

RS£C(2) 

RSEC (3) - - - 

- RSEC (NMAX) 


The remainder of the data is furnished by the flow field 
program from TAPE1 and by the mesh generation program from 
TAPE 11. These consist of: 

READ 1 

TAPE1; Q array, TAU, NC 
READ 2 

TAPE11; JMAX, KMAX, LMAX, KLMAX, JLE , JTE , KTIP, 
NBLADE , PERD, BETA34, BLDIAM 

READ 3 

TAPE11; X(KL, J) , Y(KL,J) , Z(KL,J) , KL = 1, KLMAX, 
and J = 1, JMAX. 


7.4. Program Output 

The generated output is described in this section. 

A description of the output format is provided in the next 
subsection with an example furnished in Appendix B.3. 
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7.4.1. Output Format 

To give the user a cursory look at the solution, printer 

plots of the pressure coefficient, C , are presented. These 

consist of the chordwise C variation on both the face and 

P 

camber side of the blades at varying constant radii from the 

centerline of the nacelle and of the axial C variation for 

P 

different meridional angles § on the nacelle surface . 

Next the local geometric and flow field properties of 
the face and camber side of the blade at a constant radius 
are printed. These consist of the z, <j>, and y coordinates 
at a given point and of the pressure, density, axial velocity, 
radial velocity, rotational velocity, and energy at the same 
point. 

The final major printouts consist of 

• the index of each station 

• the radius measured from the centerline of the 
nacelle 

• the chord at each section (constant radius) 

• the geometric angle at each section 

• the effective pitch angle at each section 

• the effective angle of attack at each section 

• the lift coefficient per section 

• the drag coefficient per section 

• the relative velocity per section 
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• the dynamic pressure per section 

• the thrust per blade per section 

• the torque per blade per section 

• the total thrust per section 

• the total torque per section 

• the total thrust coefficient 

• the total torque coefficient 

and of the geometric and flow field properties of the given 
blading. These consist of the number of blades, Mach numbers, 
blade diameter, free stream pressure, free stream density, 
free stream velocity, advance ratio, thrust coefficient, 
torque coefficient, power coefficient, and blade efficiency. 

7.5. Operating Procedure and Trouble Shooting Hints 

This particular program is extremely simple and straight 
forward to use. The only problems that could arise are 
possible bad data on the input cards or files . 
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8. PROGRAM SUMMARIES 

It is intended here to discuss the compatibility of the 
complex of computer codes to various computer systems. Also, 
the typical run times on a CDC 7600 machine are presented 
for time comparisons. 

8.1. Machine Compatibility 

This particular complex of codes was developed on the 
CDC 7600 machine at NASA- Ames . It was intended to create 
the codes such that they are all machine independent. A 
version of these codes are running on the IBM 370/3033 
system at Lewis Research Center. This shows the ability 
of the codes to be run on various systems . 

Typical storage requirements on a CDC 7600 are presented 
next. For the mesh generation program, which is overlayed, 
the first overlay determines the shape of grid domain and 
needs 54401 octal words of small core storage. The second 
overlay which determines the patch coefficients uses 22612 
octal words of small core storage and 360750 octal words 
of large core storage. The clustering of the grid points 
is done in the third overlay and uses 13407 octal words of 
small core storage. Finally, the fourth overlay creates the 
grid and uses 47422 octal words of small core storage and 
1017260 octal words of large core storage. 
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The 3-D flow field code uses 105737 octal words of 
small core storage and 1044340 octal words of large core 
storage . 

SCM 36650 octal words of small core storage and LCM 546460 
octal words of large core storage are used in the data 
reduction program. 

The 3-D flow field code used all available storage 
space on the NASA- Ames 7600, but if a virtual memory machine 
was available, the user could increase the useable core such 
that a far better grid resolution could be obtained. This 
should have a significant impact on the solution process. 

Even though it would take more computational time the 
obtained solution could be refined enough to give significant 
results since we believe that at present the solution 
procedure is correct, but there are not enough grid points 
available to obtain a well behaved solution. 

8.2. Runtime Estimates 

The mesh generation code and the data reduction code 
have insignificant running times as compared to the 3-D 
flow field code. Typically, these codes consume less than 
30 seconds on a CDC 7600 machine for a 45 (axial direction) 
x 21 (radial direction) x 11 (meridional direction) grid. 

There are three times given for the 3-D flow field code 
on a 45 x 21 x 11 grid system. The first time is for the 
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inviscid flow past the nacelle only, followed by the inviscid 
flow past the nacelle-blade combination and finally these 
same conditions calculated by the diagonal algorithm. 

Table 1 shows that the diagonal algorithm is 26 % faster 
than the regular algorithm, where the comparisons are based 
on solutions which were iterated from a uniform flow initial 
condition. Since only the timing ratio is of interest in 
Table 1, there was no need to run the solution to convergence. 


TABLE 1 


CONDITION 

(Inviscid) 

ALGORITHM 

ITERATIONS 

TIME 

(CDC 7600) 

RATIO 


Nacelle 

Reg. 

201 

1364 seconds 

6.79 

Nacelle- 

Blade 

Reg. 

100 

700 seconds 

7.0 

Nacelle- 

Blade 

Diag. 

300 

1545 seconds 

5.15 
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Appendix A 
Sample Input Data 

Sample input decks have been reproduced and are 
presented in this appendix for the mesh generation 
program, the three-dimensional flow field program, and 
the data reduction program. 

A.l. Mesh Generation Program 

Input data for the mesh generation program is obtained 
from the INPUT file (TAPE2) and from an archives file 
(TAPE5) containing node clustering information along 
with the geometric representation of the prop-fan 
configuration. Sample data read from INPUT is given in 
Figure A-l, while sample data read from the archives file 
is given in Figures A-2 and A-3. 

A. 2. Three-Dimensional Flow Field Program 

Sample input data for the three-dimensional flow field 
program is presented in Figure A-4 along with a set of 
control cards applicable to the CDC 7600. 

A. 3. Data Reduction Program 

Sample input data for the data reduction program is 


given in Figure A-5. 



T 

F 

t- 

T 



X 

1 ! 

0 

27 US 

11 21 

1 1 8 

van l 

0.0 


0.5 

.02025 

.05 

* 10.0 

10 ' 0 

o.o 

.05 

■ i • 0 y 

*•9 

O.o 


.0 

- 1.0 

C.6 

C.3 


.1 

- 0.5 

* .0 

0.0 


.05 

0,00 

Uo 

O.o 


.2 

• 0.1 

0.6 

C.<l 


.2 

• 0,1 

1.0 

0.0 



.0202$ <3,5 


Figure A-l. Sample Input for the Mesh Generation Program- -Read from Input 
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45 

-,?b°0 
-.267713 
-.266025 
-. 25^35 

-.251 ?439 

-,?«1A53 
-.23*46 
-.22*072 
-.221681 
-,2l"?9 
- . 2 Of 0 
-, 1 3 

- . 1 8>- 066 
•.17 C 6S 

- . 1 5 4 8 1 7 
- . 1 4 '- 4 
-.4209413 
0.10144 
.112688 
.1356 
. 1 62^92 
.1^41352 
.225 3 93 
,25o' 36 
»272' 4 
,2<ll7(,q 
,31 I *55 
,330'- 
.35 O' 4 5 
,36'MO 
. 388*35 
,4?7^?5 
,4667 1 5 
.505605 
,52505 
,563 c 4 

,6349 

,705uQ 
.7756A 
, A 4 62 A 
.91673 
,98 7 3 p 7 
1 .05* 


0.0 

00 3 

0 0 C 9 8 <j 

0 1)7076 
C 1 4 1 2 1 a 
OI9933 
025*54 
029033 
0 3 4 0 1 5 5 
0 3 7 / 4 k 5 
04 1 V 75 

0 4 4 04 7 

04 8 3133 
0 5 2 6 7 7 7 
9 5 9 1 5 0 
OoO 35 8 4 

0 9 4 » 9 8 7 
«C7 792 

1 1 1 ”04 1 
155*2 
159*35 

1 06 11 0 3 
W,o«2 
1 7 6 2 7 j 3 
1 7 0 1 5 7 

1 8 0 *’ 8 

1 8 J "44 
18 1*7 
» 80^32 

1 8ri73 
1 7q 5 i4 
1 7 8 9 4 4 
1 77 1 44 

1 7 4 u 0 9 
1 7n«02 
1 6 7 0 0 5 
1 ft c e 7 e 
1 opOPH 
1 5 A 2 4 4 
157*26 

1 57‘ , o40 
1 60 On 
1 r> p 7 IP 
1 65 ->74 
167*3 


Figure A-2. Input to Mesh Generation Program from 
Archives File- -Nodal Clustering Data 
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ti 


.0 


t 1 


S3 ]i 25 


2°. 

ft 

,o5ae 

.1395 





1 

°.o 

-0,0000 CO 

-o.oooooc 

o.c loo 

0.0l37<j9 

-0,619529 

0,2392 

? 

°,o2o<' 

0,01*60* 

-o.o2oi 3 * 

0.0300 

0 , 0222 1 9 

-0,025)61 

0.2392 

3 

o.o y oo 

0.025593 

•0. 02^5*5 

o.cSoo 

0.022751 

-0,013539 

0.2392 


n • 0 9 0 0 

p . 0 3 1 * H a 

•0.o3 7 256 

0. n 7 0O 

0. 03«393 

•0.090793 

0.2392 

5 

a.o«oo 

0.036917 

-0.0« J 163 

0.C900 

0.039287 

-o.o«7iea 

0.2392 

« 

o , i o o r 

0.0915a? 

-O.oSOOfeft 

0.1250 

0,0968 j| 

•0,057*95 

0,2392 

7 

o. iSon 

0,051660 

•0,o693qJ 

0.1750 

0 , oSe 1 06 

•0,070792 

0.2392 

P 

w, 200ft 

0 ,0603aa 

-0.0 7 68jr 

0.2250 

0.09935 1 

-0,0*270! 

0.2392 

9 

0.2^00 

0,06*051 

-0.0 fi,, 37C 

0.2750 

0, 071557 

-0,093803 

0.2392 

l n 

o.joor 

n , 0 7 a ft 1 ? 

-0.0 <5 *9o<; 

0.3250 

0.077*06 

-0,103»9S 

0,2392 

1 i 

o,3Sof 

n,0?i)5c3 

-0 . 1 0*us? 

0.3750 

0,092**2 

•0,112562 

0 . <392 

1? 

>’.ai»oc 

O.OPaft n 3 

-0.1 I605 5 

0.9250 

0,026399 

•0,116902 

0,2392 

13 

o.uSoc 

0.067U13 

•0.l209fl3 

0.9750 

0,0977*3 

•0,122236 

0,2392 

1 a 

0 . 5 0 0 0 

0,067313 

-0. j225p9 

0.5250 

0, 095903 

■0,121979 

0,2392 

IS 

O.550C 

0 , 0*37*2 

»0.l20Oj^ 

0.5750 

0,0P10 7 9 

•0,117972 

0,2392 

Ift 

O.ftOor 

o,o7 M OPft 

-0.119*97 

0.6250 

0.0750*0 

•0,110*55 

0.2392 

1 7 

o ,fc5 or 

0,072137 

-0.105873 

0,6750 

0,069382 

-0,100359 

0,2392 

IP 

o.roor 

n .Ofteftfta 

-o.o^aj | 7 

0.7250 

0,060621 

-0,0*7165 

C.2392 

19 

U,7Spr 

0,062*72 

-0.O 7, 59? 

0.7750 

0,060073 

-0,071599 

0.2 392 

?0 


0 . 0 5 7 " ft 9 

-0.0*316* 

0.7250 

o.o332o7 

•0,059*11 

0 , ? 392 

2 1 

0 . ft 5 n i 

0 ,0«5'*C6 

-0.09*10* 

0.6750 

a.092* t e 

-c,o3eo£n 

0.2392 

£? 

",<)U pi 

0 . 0 3 " * a 7 

-0,n30ftP 1 

0.9100 

0,033672 

-0,027961 

0,2392 

23 

0 .9 2 on 

0.0307efi 

-0.n253 h 3 

0.9300 

0, 027755 

-0,022fl3« 

0,2392 

2 0 

0,5am 

0.02 a 6*2 

-0.020337 

n , 9 5 n 0 

0.0219(12 

-0,017813 

0.2392 

25 

0,9ooft 

0.01*167 

-0.01522ft 

0.970O 

0.0198,3 

-0.012567 

0 , <39 2 

2ft 

0.9*00 

0,011369 

-O.o09*85 

0.9900 

0,007*39 

-0,007179 

0,c392 

27 

1,0000 

n . 0 

0.0 




0,2392 

1”. 

P 

. 0 * 3 6 * 

. 1 am 





1 

0.0 

•rwOOOOCO 

•o.ooooof 

0 . 0 1 00 

0,0137*9 

•0,019529 

0.2392 

? 

<■ . 0 2 0 (' 

o.oi*ftc* 

-O.o20 1 3* 

0 . 0 3 0 0 

0 , 0 2 2 2 1 9 

-0 , 025 1 6 1 

0,2392 

3 

0 . 0 a 0 0 

0.025593 

-0.0295ft5 

0,0500 

0.02*75! 

•0.033539 

0 . « 392 

0 

0,0*00 

o.0316*a 

-0 . o3 7 256 

0 . ft 7 0 f' 

0.0)03*3 

-0.090793 

0 . ? 3 9 £ 

5 

0 . n « 0 1 

0,0169 J7 

-0.0“916 7 

0.0900 

0,039287 

-0,0973*2 

0,2392 

ft 

I’.lOoo 

0 . oa 1 5a2 

-U.o509 6t 

0.1250 

0.0068, 5 

-0,057695 

0 , e 39 2 

7 

O . 1 So ( 

0.051660 

-0.06939? 

0.1750 

0 ,056 1 <(6 

-0,0 7 0792 

0,2392 

« 

0,2001 

n ,0fc03aa 

•0. 0^**3 0 

0.225 0 

0,0693,1 

-0,0*2708 

0,2392 

9 

0.2S00 

n. 06*051 

-0.06*370 

0.2750 

0,071557 

-0,093*03 

0.2392 

1 0 

0.300c 

0.0»aRi2 

-O.o <, *9g5 

0.3250 

0,0778 0 6 

•0,103898 

0,2192 

1 1 

o.35or 

0»0*05c3 

-0. 10*95? 

0.3750 

0.0*2*62 

•0. 112592 

0.2392 

1? 

o.aooo 

n.0*a*a3 

-0.1 16059 

0.9250 

0,096399 

■ 0,1 18902 

0.2392 

1 3 

0 . a * 00 

o.087u 1 3 

.0. 1 2098 t 

0.9750 

0.0977*3 

•0,122236 

0.2392 

1 0 

a. 5001 

0.0*7313 

-0. 1 22585 

0.5250 

0,095<)0S 

-0.121979 

0.2392 

15 

0 . 5 5 0 c, 

0,0»37P2 

-O.i20«i* 

0.5750 

o,oe 1079 

•0,117972 

0.7392 

1 ft 

, b 0 0 0 

O .07«Oft6 

-0. 1 1«697 

0.6250 

0 , 0 7 5 n * 0 

• 0,1 10655 

0.2392 

1 7 

a .ftSoo 

0.072137 

.0 . 1 05*71 

0.6750 

0,009382 

.0,100359 

0.2392 

1 P 

0.700c 

0.066*60 

-0 . 0^9 1 j 7 

0.7250 

0.069621 

•0.0*7165 

0.2392 

1« 

o,75 on 

n.06207 ? 

-0.0 7t> 59? 

0.7750 

0.060(173 

•0.071599 

0.2392 

2n 

0 , P 0 0 ( 

0.0570*9 

-C.063l6« 

0.8250 

0,o532 0 7 

•0,0*9*11 

0.2392 

?! 

0 ,65or 

O.OOftOQfi 

-0.096105 

0.6750 

0,0028,8 

-o.oieoOo 

0.2392 

2? 

0.900c 

O.03O007 

— 0.o3C6g1 

O.olOO 

0.033*72 

•0,027*61 

0.2392 

23 

0 , 9 2 0 0 

O,03" 7 6* 

- 0.0 253&7 

0.5300 

0,027755 

-0,022*38 

0.2392 

2a 

0,9a gn 

0,02*602 

-C.o2f33 7 

0 . 9 5 0 0 

0, 021002 

-0,017813 

0.2352 

25 

0,9600 

0 . 0 1 * 1 6 7 

-0.olS??f 

0.9700 

0.0198,3 

-0,012567 

0,2392 

2ft 

0.9*00 

0.011 3o9 

-O.o09*h5 

0.9900 

0,007Pj9 

• 0, 007)79 

0.2392 


Figure A- 3. Input to Mesh Generation Program from Archives 
File- -Geometric Representation of Prop-Fan 
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Figure A-3. (Continued) 
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Figure A-3. (Concluded) 
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Figure A-4. Sample Input for Three-Dimensional Flow Field Code Including 
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Appendix B 
Sample Output Data 

Sample output has been generated and is presented in 
this appendix for the mesh generation program, the three- 
dimensional flow field program, and the data reduction 
program. Typical outputs from the first two programs are 
very lengthy and are presented here in a very abbreviated 
form. 

B.l. Mesh Generation Program 

Sample outputs are presented in Figures B-l, B-2, 
and B-3 showing the Nacelle coordinates , the definition 
of one of the patches, and the coordinate locations in 
both Cartesian and cylindrical coordinates. 

B.2. Three-Dimensional Flow Field Program 

Sample outputs are presented in Figures B-4, B-5, and 
B-6 showing functions of time, surface variables, and 
general flow field variables. 

B.3. Data Reduction Program 

Sample outputs are presented in Figures B-7, B-8, 
and B-9 showing the printer plots obtained for C^, a 
summary of the local properties, and a summary of the 
force data. 
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Figure B-l. Sample Mesh Generation Output- -ITacelle Coordinates (z,R) 
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Figure B-2. Sample Mesh Generation Program Output-Patch Group Definition 
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Figure B-3. Sample Mesh Generation Program Output- -Cartesian and Cylindrical 
Definition of Each Coordinate Point 
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Sample Flow Field Output--Function of Time 
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Figure B-5. Sample Flow Field Output--Surface Variables 
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Figure B-6. Sample Output- -Variables Across the Flow Field 
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Figure B-9. Summary of Forces from the Data Reduction Program 
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